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I. INTRODUCTION 


Functional integration is an extremely useful tool for the non-perturbative analysis of 
quantum field theories (QFT from now on). To ensure that calculations involving these 
integrals be well dehned, it is convenient to regularize them by working on a spacetime 
lattice. A continuum QFT is then defined in terms of a sequence of lattice held theories. 
There may be many lattice representations for a single continuum theory, so it is not obvious 
that such a dehnition is unique. In the majority of lattice held theory simulations, the issue 
of uniqueness is addressed by appealing to the principle of universality. Universality is the 
property that, in the continuum limit, the system ceases to depend on the exact nature 
of the lattice theory, and many diherent lattice held theories lead to the same QFT. In 
the approach to the continuum limit, which is located at special—critical—values of the 
coupling parameters, various quantities are expected to have a universal dependence on the 
couplings. Two diherent lattice theories that display the same behavior at this critical point 
are said to belong to the same universality class. 

In the present paper, lattice models that do not to conform to this expectation of uni¬ 
versality are described, and the extent of universal behavior is clarihed somewhat. Four 
models are considered here: two lattice versions of the S'0(1,1) nonlinear sigma model and 
two lattice versions of the SO{2, 1) sigma model. In each case the two lattice models become 
the same classical theory in the classical continuum limit. 

Numerical evidence strongly suggests [^-|^ that both lattice 0(2,1) models do belong 
to the same universality class in the sense that they both describe a pair of free uncoupled 
massless scalar helds in the continuum limit (they are trivial). The situation for the 50(1,1) 
models is not as clear. One lattice 50(1,1) model is known to become a free massless scalar 
held in the continuum limit, but the continuum limit of the other lattice model is unknown. 
The current numerical data does not rule out triviality. 

All of these lattice theories, however, have an order parameter, called the magnetization, 
which is sensitive to the details of the lattice model near the critical point. This is unlike the 
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behavior observed in more familiar models, such as the Ising model or the compact nonlinear 
sigma models IQ. It is argued that the cause for this nonuniversal behavior is the radically 
different scaling of the helds with a parameter (3 (analogous to an inverse temperature in 
statistical mechanics) that is made possible by the noncompact conhguration spaces. This 
point will be further discussed in Section 0. 

The remainder of the paper is organized as follows. In Sec. |1|, the nature of the approxi¬ 
mation involved in dehning a QFT on a lattice is explained, and the freedom involved in this 
approximation is demonstrated. The relevance of universality to lattice held theory is then 
explained. The language of the renormalization group (RG from now on) is used to show 
why certain quantities {e.g. critical exponents) are expected to be the same for lattice held 
theories that are in the same universality class. In addition, the process through which two 
diherent lattice theories become the same QFT in the continuum limit is explained in terms 
of the concept of irrelevant operators. The general structure of nonlinear sigma models is 
given in Sec. Sec. ^ introduces two possible lattice 5'0(1,1) models and shows that 
they have diherent behavior in the critical region. Sec. |V| discusses the analogous situation 
for the SO{2, 1) model. In Sec. [V| the reason for the nonuniversal behavior is discussed, 
and the conclusions are presented. 

It should be noted that unlike the two and three dimensional models that are of interest in 
the related discipline of critical phenomena, this article deals with models in four spacetime 
dimensions, since the motivation is to study quantum held theory. To carry out numerical 
simulations of the lattice held theory, it is necessary to analytically continue the time into 
the Euclidean regime. The result of this Euclideanization is mathematically equivalent to a 
statistical mechanical partition function, with a Hamiltonian given hj -^ = Se where 5^ is 
the Euclidean action of the lattice held theory. 
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II. LATTICE FIELD THEORIES AND UNIVERSALITY 


A quantum field theory is constructed from a lattice field theory as follows. The QFT 
is dehned by the vacuum expectation values of time ordered products of operators built 
out of the local quantum helds. These time ordered products can be expressed in terms of 
functional integrals. The functional integrals are in turn dehned by discretizing spacetime, 
which turns them into multiple integrals, albeit in a very high dimensional space. The 
multiple integrals are then amenable to numerical analysis through Monte Carlo techniques. 
The lattice held theory approximants to the expectation values of the continuum QFT can 
then be calculated^. The simulations are carried out on successively larger lattices, and the 
approximants are then extrapolated to the continuum limit to yield the desired continuum 
vacuum expectation values. 

Thus, there are two identihable stages in the process of quantization through functional 
integrals. First, a lattice held theory is dehned as a large dimensional multiple integral with 
an action constructed from the classical action. Second, the continuum QFT is recovered 
as the limit of a sequence of lattice theories simulated on successively larger lattices, with 
the lattice spacing simultaneously considered to be taken to zero. If the continuum QFT 
contains any hnite length scale, such as the inverse of the mass or the correlation length 
this length scale becomes inhnite when measured in units of the lattice spacing a, which 
goes to zero in the continuum limit. Therefore, in the second stage, when the lattice spacing 
is taken to zero, it is necessary to identify a point in the parameter space of the theory 
where the lattice version of the characteristic length scale of the QFT diverges in units of 
the lattice spacing. The correlation length of the nonlinear sigma models discussed here is 
already inhnite on an inhnite lattice because the models are massless. Instead, the relevant 
length scale here is which is associated with the interaction. The continuum limit is 
sought where the characteristic length scale measured in units of the lattice spacing. 




^Of course, non-numerical techniques may also be used. 
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diverges {j3R = This means that a continuum limit exists only where (3^ 0. 

The lattice action constructed during the first stage is not unique, because the continuum 
action contains derivatives of the helds. The customary prescription for passing from the 
classical continuum action to the lattice action is to replace i9^0(a;) — ^Afj_(l){x), where 
A^0(a;) = 0(a; + /r) — 0(a;). Here x = {ai, aj, ak, al) where a is the lattice spacing, i, j, k, I G 
,N} are integers, and a; + /i (/i G ,4}) denotes the lattice site which is the 

nearest neighbor of x in the forward /i direction. In the classical theory, one can perform a 
point transformation 4>{x) <h(0(a;)), and the equations of motion in terms of the new helds 

<h(a;) have the same content as the equations of motion in terms of the old helds. However, 
the lattice action obtained from the classical action expressed in terms of the helds <h(a;) will 
be diherent from the lattice action obtained from the classical action expressed in terms of 
the helds 4>{x). The equivalence of the theories at the continuum classical level rests on the 
chain rule for diherentiation, which is not valid on the lattice. It is not a priori clear that 
the continuum QFT obtained at the second stage is the same for the two diherent lattice 
actions. 

In the great body of lattice held theory simulations, this possible ambiguity is circum¬ 
vented by appealing to the principle of universality. The principle of universality maintains 
that, in the critical region, the large distance behavior of lattice systems is mostly inde¬ 
pendent of the small distance features, such as the detailed choice of the lattice action. In 
extrapolating to the continuum QFT, one reexpresses the observables of interest in units of 
appropriate physical lengths. In the continuum limit, as noted above, these hnite physical 
units of length become inhnite as measured in terms of the lattice spacing. Therefore, the 
specihc features of the lattice action automatically become small distance features as com¬ 
pared to anything of interest in the continuum QFT. Universality requires that these small 
distance features be unimportant. 

More precisely, all possible lattice models separate into classes. All lattice held theories 
within a single universality class become the same QFT in the continuum limit, but lattice 
models belonging to different universality classes yield different continuum QFT’s. From 
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the point of view of quantum field theory, the details of the lattice theory do affect a small 
set of quantities that are nonuniversal in the critical region, but which do not affect the 
physical content of the continuum QFT. On the other hand, in statistical mechanics, these 
nonuniversal quantities may have physical content. For example, in the models considered 
in this article, it is found that the magnetization, which is related to an expectation value of 
a field (Eq. [4.41) in one of the formulations of the SO{2, 1) model, has nonuniversal behavior. 
The magnetization is relevant for the wavefunction renormalization [Q], but apparently does 
not affect the S matrix of the theory (which is trivial). 

There are three criteria which determine a universality class. They are ( i P- 80): 

1. The spacetime dimension of the system; 

2. The internal symmetry group of the system, and which representation of this symmetry 
group is furnished by the helds in the lattice action; 

3. The nature of the critical point. 

These three criteria alone still leave one with the freedom to use many different lattice actions 
to simulate a given theory. Although the above criteria are widely known, additional criteria 
may be needed to determine a universality class uniquely, cf @ p. 18. 

It is now seen how universality leads to a resolution of the ambiguity inherent in formu¬ 
lating a lattice held theory. To ensure that the lattice held theories under consideration go 
to the same continuum limit, the lattice actions must have the same spacetime dimension, 
have the same internal symmetry group, and be constructed out of helds belonging to the 
same group representation. It is not immediately clear that diherent lattice actions (such as 
those obtained from the helds 0(x) and d>(x) in the discussion above), which satisfy criteria 
1) and 2) will satisfy the third as well. It is commonly assumed that the hrst two criteria 
are sufficient. This is occasionally tested when numerical studies are carried out. This as¬ 
sumption is often made for scalar helds in four spacetime dimensions, and is buttressed by 
the strong evidence existing [[^] for the triviality of A<h^ theory. However, for other scalar 
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field theories, as yet unknown fixed points may exist. 


The principle of universality appears most natural in the framework of the renormaliza¬ 
tion group (RG). A good introduction to the subject may be found in P, g], [Q or 
while a detailed review is given in and [Q . The following is a brief overview of the key 
ideas. As noted previously, the continuum limit is extracted from the features of the lattice 
theory with length scales much larger than the lattice spacing. The RG provides a means 
of studying these features by gradually integrating the high momentum degrees of freedom 
out of the lattice version of the partition function. 

A renormalization group transformation consists of the following steps: 

1. Begin with an action dehned at the scale of some cutoff A (on the lattice A ~ i). This 
action should be the most general possible action consistent with the symmetries of the 
theory] 


2. Integrate modes with momenta in the range ^ < k < A, out of the functional integral. 
Here s is some scale factor greater than one; 


3. Rescale all the remaining momenta by fc —>■ sfc so as to restore the original cutoff A; 

4. Redehne the dummy held variables of integration in the partition function by cf'ik) = 
a{s)(j){k), where a{s) is the function required to restore the kinetic term in the action 

to its original normalization. 


The effective action obtained after an RG transformation has the same form as the original 
action, but the coupling parameters are modihed. RG transformations thus generate a how 
in the space of coupling parameters. In the cases of interest, this how terminates at a hxed 
point, at which the ehective action is RG invariant. 

If the parameters in the original action have their critical values and the correlation 
length ^ is inhnite, then ^ remains inhnite under RG transformations. The hows that leave 
^ inhnite generate a subspace of the coupling parameter space known as the critical surface. 
It is assumed that every point on the critical surface hows to some hxed point. The theory 
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at this fixed point has the same long distance behavior as any critical theory connected to it 
by RG transformations. It should be emphasized that the RG transformation changes the 
natural unit of length (the lattice spacing). The long distance behavior is only the same in 
both cases when the same unit of length is used. Therefore, the long distance behavior of 
the original critical theory can be deduced by studying the theory at the fixed point. 

The origin of universality now becomes clear. The critical surface breaks up into domains, 
such that all the points in a domain flow to the same fixed point. Each domain defines a 
universality class. Within a given universality class, the details of the lattice action of some 
starting critical theory are not important because all theories have the same long distance 
behavior as the fixed point theory. 

As will be shown below, the long distance behavior for a theory which is slightly off the 
critical surface, is the same as the behavior of the theory at some point in the neighborhood 
of the fixed point. For simplicity, the analysis of RG transformations in the neighborhood of 
the fixed point will be restriced to the case of a single coupling parameter. Under a general 
RG transformation, 

K' = R{s,K) (2.1) 

the function R{s, K) is a smooth function of the coupling parameter K for a finite rescaling 
of the cutoff by s. If K is sufficiently close to the fixed point, where the coupling parameter 
has its fixed point value K*, the general RG transformation is well approximated by the 
linear part of its Taylor expansion: 

K' = K* + {K - K*) = K* + k{s){K - K*). 

K=K* 

Gonsidering two successive RG transformations with scale parameters si and S 2 , one finds 
K" = K* + A{si)A{s 2){K - K*) = K* + A{siS2){K - K *). 

This is just the semi-group property of the RG. The only possible functional form for A(s) 

is 
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A(s) = . 


( 2 . 2 ) 


It is now seen that the behavior of the coupling constant K near the hxed point under 
RG transformations is determined by the exponent y. \i y > 0, the coupling parameter 
is enhanced under an RG transformation, and K is called a relevant coupling. If ?/ < 0, 
it is suppressed, and it is called an irrelevant coupling. y = 0, the coupling remains 
constant under RG transformations (at least in the linear approximation) and it is then 
called marginal. In the general case of an infinite set of couplings Ki, A(s) becomes an 
inhnite-dimensional matrix. The coupling parameters no longer transform as simply as 
in Eq. (|2.2| ). Instead, the eigenvalues of A(s) transform in this simple manner. These 
eigenvalues are usually called scaling fields. The coupling parameters, however, are linear 
combinations of the scaling fields. The operators multiplying the scaling fields in the action 
are also called relevant, irrelevant, or marginal, according to the behavior of the associated 
scaling helds. Terms multiplying irrelevant scaling fields in the effective action become small 
as successive RG transformations drive the system to the fixed point. 

If the system is initially slightly off the critical surface, the direction of the RG flow of 
the relevant scaling fields will be away from the fixed point (and hence the critical surface), 
while the direction of flow of the irrelevant scaling fields will be towards the fixed point. The 
relevant scaling fields must vanish on the critical surface. Otherwise, RG transformations 
will take the system out of the critical surface. Therefore, from a starting point sufficiently 
close to the critical surface, RG transformations first drive the system to the neighborhood 
of the fixed point, because the irrelevant operators dominate. The long distance behavior 
of the original critical region is thus seen to be the same as the long distance behavior in a 
neighborhood of the fixed point, since flow trajectories connect the two regions. Eventually, 
the relevant scaling fields will drive the system away from the neighborhood of the hxed 
point. 

The RG also explains why the critical exponents of certain quantities are universal. 
Their critical behavior can be analyzed in terms of the behavior near the hxed point. Such 
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quantities themselves obey transformation laws similar to those of the coupling parameters: 

A' = sy^A. (2.3) 

One such quantity is the correlation length, whose scaling behavior is 

e(iri, 0,..., Li, L2, ...) ^ eiK 0,..., L;, L',...) = 0,..., ...), 

where Ki are the relevant scaling helds and Li are the irrelevant scaling helds. The system is 

considered initially to have all the scaling fields except Ki set to zero, and 1, since the 

1 

initial point is very close to the critical surface. Performing a transformation with s = , 

one gets 

1 _£1 _£2 
L[) = 0,..., iPi Li, iPi L2, ...) • 

1 

Since 1, and x* < 0 for all i, it follows that 

L[) ^ 0, 0,..., 0, 0,...) ~ K;^ . 

This means that the critical exponent for the correlation length with respect to the coupling 
Ki is The critical exponents of ^ with respect to any other relevant coupling param¬ 
eter can be obtained by starting with that coupling as the only non-zero relevant coupling 
parameter. The reason that the critical exponents are universal now emerges. The values 
of the critical exponents are generically determined by the eigenvalues of the linearized RG 
near the hxed point in the given universality class, and not by the values of the coupling 
parameters themselves. Any system in a given universality class will therefore have the 
same critical exponents with respect to given relevant coupling parameters (such as the 
temperature, the external held, etc.). 

III. NONLINEAR SIGMA MODELS 

Numerical results for the nonlinear sigma models with noncompact conhguration spaces 
demonstrate the need for a closer look at universality. Nonlinear sigma models are a class 
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of scalar field theories that can be obtained by imposing a constraint among a set of free 
massless scalar helds 0. The models are constrncted to be invariant nnder the global 
transformations of some hnite-dimensional Lie gronp, G. The conhgnration space in which 
the helds live is generally a coset space of G. For the compact 50(2) and 50(3) sigma 
models, the conhgnration spaces are the circle and the 2-sphere respectively. These sigma 
models are analogous to spins systems in statistical mechanics. If G is a noncompact group, 
then the conhgnration space may also be noncompact. This article will be concerned only 
with Euclidean models, invariant under SO{n,m) groups. 

Consider a set of /c -I- 1 = m + n massless free scalar helds whose Euclidean action is 
invariant under transformations of G. 

S[^p] = ^ J d'^xpij i,j e {l,2,...,k + l}. (3.1) 

The helds are cartesian coordinates for the internal space associated with each space-time 
point. In this coordinate system, rjij is diag(l, 1,...,—1,—1,...). 

n m 

If the helds at every point are required to satisfy the constraint, 

- Vij (3.2) 

then the theory described by (|3.1|) together with (|3.2|) is known as a nonlinear sigma model. 
The constraint merely requires that the length of the held vector (in the metric of the 
conhgnration space) be a constant. Since both 5[(p] and the constraint are invariant under 
the global action of G, the sigma model clearly possesses the same symmetry. Diherent 
groups describe diherent theories, and the symmetry is generally used to distinguish among 
them. 

The helds are constrained to a /c-dimensional curved surface, G, in the (/c-|-l)-dimensional 
hat space. This surface is a space of constant curvature and is generally the manifold of 
a coset space of G. The space G may be noncompact, if G is a noncompact group. The 
curvature of the conhgnration space introduces interactions into the theory. The nature 
of these interactions is most easily seen if the constraint is used to eliminate one of the 
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fields from the theory. After eliminating the action can be expressed in terms of the 

remaining helds and Gabi’^), the induced metric on the constraint surface. 

1 


SM = j d*x a, be {1,2,..., A^}, 


(3.3) 






’ 

Expanding Gab{^) around (p°‘ = 0 gives 


a, b, c, d,r,s E {1,2,..., k}. 


S[(p] = 2 7^^^ - II ^ 


VcLcVbdVrs^ 9 ^ 9 ^ ^^ 9 ^ ^^ 9 ^ 


c, ^d,^r ,^s ‘ 


(3.4) 


It is clear from (|3.4|) that plays the role of a coupling constant. It is also seen that 
the interactions involve held derivatives and that the action is not polynomial in nature. 
In four dimensions, /i has dimensions of mass and the coupling constant has dimensions 
of [mass]“^. Well known power counting arguments indicate that theories with couplings 


having negative mass dimension are not perturbatively renormalizable |^. In addition, 
naive scaling arguments indicate that that l//i^ will be irrelevant under the renormalization 
group transformation described in Sec. |I|. For this reason, it is generally believed that the 
nonlinear sigma models in four dimensions are all in the same universality class as the free 
held theory (they are trivial). 

For purposes of numerical simulation, it is convenient to render the helds dimensionless 
by the rescaling 0“ = j[i. In terms of these dimensionless helds, the sigma model action 
becomes 


= y / d'l a,b e [1,2,... ,k}, 


(3.5) 


where 


CAd 


Gab{4>) = Vab - 


Vacriad (t)''(t>' 
1 + hr- 


a, b, c, d, r, s G {1, 2,..., k}. 


When spacetime is taken to have periodic boundary conditions, any constant held is 
a solution of the classical dynamical equations which minimizes the Hamiltonian. The 
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existence of a continuous global symmetry dictates that sigma models have a k-parameter 
family of ground state solutions. Because of this, the models display the phenomenon of 
spontaneous symmetry breaking, in which only one of the possibilities is realized as the 
physical vacuum. 

Euclideanized nonlinear sigma models with compact symmetry groups have been widely 
studied as statistical mechanical spin systems, and the universality classes of these models 
appear to be distinguishable by the standard rules 0. As will be shown in this article, the 
situation is different for sigma models with noncompact symmetry groups. The 5'(9(1,1) 
and 50(2,1) models are the simplest sigma models of this type. Both of these exhibit 
behavior that would normally be considered nonuniversal. The 50(1,1) model is particularly 
illustrative, since it is classically related to the free theory by a transformation of held 
variables. Different discretizations, however, lead to different lattice theories, only one of 
which is simply related to the continuum free theory. 


IV. THE 50(1,1) MODEL 

The 50(1,1) Sigma Model is dehned by (|3.3|) with {k + 1) = 2 and pij = diag(l, —1). In 
terms of dimensionless helds, the action is 



(4.1) 


subject to the constraint that (0^)^ — = 1- This constraint surface is disconnected, so 

the theory will be restricted to the connected component dehned by (ff^ > 0. The internal 


space shown in Fig. |l] is the coset space 50(1,1)/(Z2 x Z 2 ). It is topologically equivalent 
to R and is intrinsically hat. Group transformations leave the surface invariant. 
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FIGURES 



The 50(1,1) Sigma Model defined in this way is classically equivalent to the theory of 
a massless free scalar field. Let (f)^ = sinh s and = cosh s in (|4.1D . The action in terms of 
the field s is 

2 

S = Y J d'^xd^sd^s, (4.2) 

which is just the massless free theory. The variable s measures the arc length from the 
axis along the hyperbola in Fig. |l[ 

For numerical simulations of the quantum theory, the 50(1,1) Sigma Model must be 
defined on a four-dimensional hypercubical lattice with periodic boundary conditions. The 
lattice consists of points (or sites), N along each spacetime direction. The sites may be 
connected by imaginary links, which join each site to its nearest neighbors. The distance 
between lattice sites is denoted by a, and the length of one side of the lattice is L = Na. 
Lattice approximations to the continuum action are generally obtained, as noted in Sec.||, 
by replacing derivatives with lattice differences and integrals by sums. When this procedure 
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is performed on the 5'0(1,1) model, the action takes the form 

■S = fi;-CW' (4-3) 

Here (3 is /i^a^—a dimensionless lattice version of the coupling. 

The exact form of depends on the details of the differencing procedure (two of 
which will be discussed shortly), but it generally contains interactions which tend to align 
the helds at neighboring sites. In this respect, the lattice sigma models are analogous to 
ferromagnetic spin systems in statistical mechanics. The held values at each site play the 
role of magnetic spins, and (3 plays the role of an inverse temperature. This analogy is more 
easily visualized for the compact sigma models, but it is useful for the noncompact models 
as well. 

If the direction of symmetry breaking is chosen so that the average value of 0^ is zero 
(this must be enforced by hand in hnite-lattice numerical simulations), then the average 
value of is a natural order parameter for this model. When (3 is large, {(jP‘) is near unity. 
This indicates that the helds at most sites are nearly aligned (highly ordered) around s = 0. 
For small (3, (0^) is large while (0^) remains zero. This indicates that the helds huctuate 
widely over the conhguration space and that the system is highly disordered. 

Evidently, (0^) is analogous to the magnetization in ferromagnetic systems, but because 
the conhguration space is not compact, it behaves quite diherently. In order to dehne a 
magnetization that behaves in a familiar fashion, the quantity 

= ( 4 . 4 ) 

will be designated the “magnetization” for this model. 

Ferromagnetic systems generally have some characteristic temperature, called the crit¬ 
ical temperature, below which the system becomes spontaneously aligned. Sigma models 
generally have a characteristic value of (3, denoted by /5c, above which the helds begin to 
become aligned. This is called the critical value of (3, or the critical point of the theory. The 
value of Pc depends on the lattice model and the dimension of spacetime, but for compact 
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models, it is generally nonzero. For noncompact sigma models, however, /3c is always zero. 
This is because the fields are always aligned to some degree unless an infinite “tempera¬ 
ture” (/3 = 0) induces fiuctuations over the entirety of the noncompact configuration space. 
For a finite temperature (finite j3), fiuctuations will be finite and the magnetization will be 
nonzero. A plot of the magnetization as a function of (3 for two lattice *S'(9(1,1) models in 
Fig. 1^ illustrates this. 

The normal lattice approximation to the integrand of (|4.2D is simply 

1 


{^^,s) 


:(A^s), Af,s = s{x + p) - s{x), 


(4.5) 


where the 4-vectors p are given by 


1 = (a, 0,0,0), 2 = (0,a,0,0), 

3 = (0,0,a, 0), 4= (0,0,0, a). 


Application of the same prescription to the integrand of (|4.1|) gives 


= a“^{[sinh(s(a; -1- p)) — sinh(s(a;))]^ 

— [cosh(s(a; -|- p)) — cosh(s(a;))]^} 
= 4a“^ sinh^ ^ 


= a 


-2 


-1 


(4.6) 


As a —>■ 0, the expressions ([4.5|) and ([4.6|) clearly become the same when the fields are 
smooth. In four dimensions, however, the fields are not smooth. One of the odd properties of 
functional integrals is that the field configurations that contribute significantly to the integral 
are generally neither continuous nor differentiable (except in one dimension, where they are 
continuous). This issue is discussed at length in |^, but qualitatively correct behavior can 
be surmised from the rule of thumb that configurations which contribute significantly to the 
functional integral are those for which, (S) ~ where d is the dimension of spacetime. 
Applying this rule of thumb to the lattice action implied by (^4.5|) gives 
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The same rule applied to ( |4.6|) gives 


A^s ~ a 


l-d/2 


A^s ~ ln(a) 


as a —> 0. This causes the lattice theories dehned by ( [1.5| ) and (|4.6|) to behave very differently, 
despite the fact that they were obtained from the same classical action. 

The numerical results shown in Fig. ^ clearly confirm that ([4.5|) and (^4.61) produce 
significantly different results for the critical exponent of the magnetization. These numerical 
simulations were performed on a lattice of 8“^ sites; at this point, accurate error estimates 
are unavailable. 



P 


FIG. 2. Magnetizations for the 50(1,1) Sigma Model. 


In fact, the (3 dependence of the magnetization corresponding to (|3) can be calculated 


analytically and is found ^ to behave as exp(—which does not allow a definition 
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of a critical exponent at all. The numerical results for (^4.6|) , however, indicate that the 
magnetization behaves approximately as (3^'^ near the critical point. There is a signihcant 
qualitative difference between these two lattice models. 


The radical difference in the behavior of the magnetizations for these two theories raises 
the possibility that the two models may belong to different universality classes; the model 
dehned by (|4.6|) may not become a free continuum held theory. The numerical simulations 
performed to date 0, do not settle this question dehnitively. The same simulations, how¬ 
ever, do show some differences in the critical exponents for the specihc heat, cl, of the 
two models (this result is further discussed in P]). For the model with the action ( [4.6D , 
Cl ~ while for (^4.5|) , cl ~ This difference, while not so drastic as the one for 

the magnetizations, appears to be another example of nonuniversal behavior. 


V. THE SO(2,l) MODEL 

The S'0(2,1) Sigma Model is dehned by (|0|) with (fc -|- 1) = 3 and r]ij = diag(l, 1, —1). 
In terms of dimensionless held variables, 

S' = y y -f , (5.1) 

with the constraint that (0^)^ — ~ i.'YY — 1- Again, the connected sheet 0^ > 0 is 

chosen. The resulting conhguration space, shown in Fig. |^, is SO(2,1)/(S'0(2) x Z 2 ). This 
is topologically equivalent to but is not hat. 

Substitution of hnite diherences into (|5.1|) produces one possible lattice action for the 
S'0(2,1) model. Another can be obtained by means of geodesic diherencing. Let I be the 
geodesic length between two held values on the hyperboloid shown in Fig. ^ Geodesic 
diherencing prescribes]^ that the lattice Lagrangian for the SO{2, 1) model be P/a?. (Note 
that geodesic diherencing of the S'0(1,1) model yields (4.6).) 


^Reasons for using this prescription are discussed in 


0 . 
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FIG. 3. Configuration Space for the 50(2,1) Sigma Model. 


If s is the geodesic length along the internal space from the point (0, 0,1) to the point 
and 9 is the angle measnred from the 0^ axis, then the geodesic lattice action 
for the 50(2,1) model is P 


f Ejcosh 

X^fl \ 


-1 


cosh[s(a: + /x) — s(a;)] cos 


2 I 9{x + /i) — 9{x) 


+ cosh[s(x + /i) + s(a;)] sin^ 


'9{x + /x) — 9{x) 


(5.2) 


Again, /3 = /x^a^ may be interpreted as the reciprocal of a bare lattice conpling, analogous 
to an inverse temperature in a statistical mechanical model. In a fashion similar to the 
50(1,1) model, the magnetization is dehned to be (0^)“^ (with (0^) = (0^) = 0). 

Plots of the magnetizations for the two lattice models described above are shown in Fig. 
For the geodesic model, the (3 dependence of the magnetization near the critical point 
is analytically soluble P] and has exponential behavior similar to that displayed by the 
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geodesically differenced 5'0(1,1) model. The model obtained by differencing (|5.1|) , however, 
has a magnetization that displays the familiar power-law behavior near /? = 0. For a lattice 
of 10^ sites, the magnetization behaves approximately as (error estimates for this value 
are presently unavailable). Again, there is a striking qualitative difference between the two 
lattice models. 
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FIG. 4. Log Plot of Magnetizations for the SO{2, 1) Sigma Model. 

Extensive numerical simulations of the geodesically-differenced 50(2,1) model indicate 
clearly that this theory is trivial in the continuum limit ^0. Somewhat less thorough, but 
still convincing numerical studies for the naively-differenced theory [Q] suggest the that this 
theory is also trivial. In this respect, expectations from the ideas of universality appear to 
be fulfilled. Nevertheless, the difference between magnetizations of the two lattice models 
remains to be explained. If there were a noncompact statistical mechanical models the 
difference would certainly be physically significant. 
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VI. CONCLUSIONS 


The noncompact nonlinear sigma models discussed here clearly display behavior that 
would ordinarily be considered a violation of universality. One possible explanation for this 
behavior is that the different lattice theories are associated with different universality classes. 
This explanation, however, is difficult to reconcile with the renormalization group arguments 
about irrelevant operators, which suggest that the different lattice theories should have the 
same continuum limit. It was mentioned in Secs. 0 and |V| that the numerical evidence 
supports (to varying degrees) the view that, as quantum field theories, these noncompact 
models are free (trivial). This view is tenuous for the 0(1,1) model (the data is also 
consistent with a nontrivial continuum limit), but is well supported by numerical results 
for the 0(2,1) theory. If it is true, it means that, under the action of the renormalization 
group, the different lattice theories flow to the Gaussian fixed point. 

In such a case, it must be that the critical exponent for the magnetization is one of 
the quantities in the theory that does not behave in a universal fashion. This opposes 
a large body of evidence, but much of that has been obtained in the context of models 
with compact configuration spaces. To understand this nonuniversal behavior, examine the 
standard argument why the magnetization should have universal power law behavior near 
the critical point: 

1. The critical lattice theory is connected by renormalization group transformations to 
another theory, which is at the fixed point. The theory at the fixed point has the same 
long-distance behavior as the original theory. 

2. Near a fixed point of the renormalization group transformation, the coupling constant 
/3 will scale by some power of s, /3 —>• /3' = s^/3. 

3. The field will also scale by s to some power, (j) —^ Since the magnetization is 

generally just {(/>), it picks up an overall factor of 
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4. 


The above behavior leads to the relationship M'{j3') = s This relationship 

implies M{f3) 13 y . 


It is the third of the above statements that is false when applied to the magnetization of 
the lattice free held and the geodesic 0(2,1) model. In those cases, the magnetization is not 
jnst the average valne of the held; it does not scale in the same fashion as the held in the 
kinetic term of the action. This means that the factor of s~'^^ in front of the magnetization 
is not correct in this case, and the argnment for a nniversal exponent fails. This shonld not 
ahect the nniversality of other critical exponents in the theory, except perhaps those that 
are related to the exponent for the magnetization. 

The fact that the conhgnration spaces for these models are not compact plays a major role 


in the observed diherences in critical behavior. In the action ([4.6|) , the hyperbolic sine differs 
greatly from A^s over the range of field valnes that contribnte significantly to the integral 
(recall that the fields which contribnte to the fnnctional integral are generally discontinnous 
||T^). This canses the helds to spread over the conhgnration space at an exponential rate 
as /d —*• 0, resnlting in an average held which scales exponentially with (3, instead of with 
the nsnal power law. This type of scaling would be unlikely if the conhgnration space were 
compact, simply because a compact conhgnration space places a limit on the size of the 
helds. 

It may seem unusual that the local scaling of the helds can ahect behavior that is normally 
considered to be due to long range correlations. While it is the long range correlations that 
cause the helds to align, the value of the order parameter is inhuenced by the local scaling 
(with j3) properties of the held. For the sigma models discussed above, the local behavior 
is sufficiently diherent to cause qualitative diherences in the behavior of the magnetization 
near the critical point. 

In this sense, the behavior of the magnetization is like that of the critical coupling (or 
critical temperature). The value of the critical coupling is known to depend strongly on the 
local details of the lattice model; it is not a universal quantity. The conclusion of this article 
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is that, for some noncompact models, the same is true of the magnetization. In held theory, 
this magnetization is not an observable quantity, and its nonuniversal behavior should not 
affect the continuum theory. As noted in Sec. however, it is signihcant from the statistical 
mechanical viewpoint. 

Regarding conclusions about the triviality of these models, there remains uncertainty 
about the continuum limit of the model dehned by (|4.6|) . More complete numerical simula¬ 
tions of this model are necessary before its continuum limit can be cited with conviction. In 
addition, it would be helpful to obtain an analytical approximation to this theory for small 
values of /?, similar to the one that supports the numerical data for the 0(2, 1) model 0. 
These tasks remain for future research. 
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